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Abstract 

We argue that the relativistic correction 6Tlc.m,. to the center-of-mass vector can 
CsJ ' lead to the approximate equality of the proton and neutron electric polarizabilities in 

' the quark model. The explicit form of 5Hc.m.. depends only on the non-relativistic 

. potential between quarks. In particular, this correction is the same for the potential 

5^ ' generated by Lorentz- vector and -scalar interactions. 

< 

^ ■ 1 Introduction 

^ ■ The electric polarizability a and magnetic polarizability /3 are fundamental structure con- 
O . stants of the nucleon having a direct relation to the internal dynamics of the particle. Their 
magnitudes depend not only on the quantum numbers of the constituents, but also on the 
i properties of the interaction between these constituents. The prediction of the electromag- 
^ I netic polarizabilities and the comparison with experimental data may serve as a sensitive 
O ' tool for tests of hadron models. 



X 



The quantities a and f3 can be obtained from the amplitude for low-energy Compton 
scattering off the nucleon. In the lab frame, this amplitude up to 0(co'^) terms reads 0] 



T = Tsorn + aujiuj2ei ■ el + /3(ki x ei) ■ (k2 x e^) , (1) 

where uji, kj, and are the energy, momentum, and polarization vector of the incoming 
{i = 1) and outgoing {i = 2) photons {h = c = 1). The contribution Tbo™ corresponds to 
the amplitude for Compton scattering off a point-like particle with the spin, mass, charge, 
and magnetic moment of the nucleon. 

The result of the most recent and most precise experimental investigation P, ^ , obtained 
by Compton scattering off the proton below pion threshold, is shown in the first line of 
Table [^. For the neutron the same method is extremely difficult : (i) Quasi- free Compton 
scattering from neutrons bound in the deuteron has to be carried out leading to sizable 
corrections at low energies, (ii) The Born amplitude Tbo™ is very small for the neutron. 
Therefore, the corresponding interference term is also small and, hence, cannot be used as 
is done in case of the proton, (iii) The low-energy differential cross section is less sensitive 
to the electric polarizability in case of the neutron as compared to the proton. 

In order to avoid these difficulties, the differential cross section for electromagnetic scat- 
tering of neutrons in the Coulomb field of heavy nuclei |^ 

%1 = ^MHZefRea {.„si„ l-f§(l-M (2) 
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has been investigated in a series of narrow-beam neutron transmission experiments. In 
p is the neutron momentum, —a the nuclear amphtude, k„ the neutron anomalous magnetic 
moment, M the neutron mass, and e the proton charge. Of these experiments only the 
Oak Ridge experiment |0] has a reasonable precision. This result is cited in the second line 
of Table ^ where the necessary correction for the second (Schwinger) term in the braces of 
(Q) has been carried out. This Schwinger-term correction was disregarded in the original 
evaluation of the experiment ||^ and reads +e^/€^/4M^ = +0.62 ||^. 

In it was shown that it is also possible to use experimental data on quasi-free Compton 
scattering on neutrons bound in deuterons at energies between vr meson threshold and the 
A peak to extract the electric polarizability of the neutron. This method has successfully 
been tested for the proton and applied to the neutron |10] , leading to the result listed 
in the third line of Table |I|. 



Table 1: Experimental electric polarizabilities for proton and neutron 



ftp 


12.2 ± 0.3(stat) ± 0.4(syst) ± 0.3(model) 


Compton H, 1 


an 


12.6 ± 1.5(stat) ± 2.0(syst) 


Coulomb 


an 


12.5 ± 1.8(stat)+J:^(syst) ± l.l(model) 


Compton H 



From Table [l| we obtain the following result for the difference between the electric polariz- 
abilities of proton and neutron 

ap — an = —0.3 ± 1.8 (3) 

which apparently is compatible with zero. 

The electromagnetic polarizabilities of hadrons have been calculated in many different 
models. Though much effort has been devoted to these model calculations, all of them can 
not be considered as completely satisfactory. In particular, there is a problem in explaining 
within a non-relativistic quark model that the electric polarizabilities of proton and neutron 
are equal to each other as suggested in (|]) . It was derived many years ago [0, [Tl|, |12[ that d 
can be represented as a sum 

2^|(fc|D|0)|2 ^ ^ 

where D is the internal electric dipole operator, |0) and \k) are the ground and excited 
states in terms of internal coordinates, Ek and Eq are the corresponding energies. Note that 
P|0) = P|A;) = 0, where P is the operator of total momentum of the particle. The term Aa 
in a has a relativistic origin and its leading term is equal to 

Aa = ^ , 5 

3M ' ^ ^ 

where Z and M are the particle electric charge number and mass, and r^; is the electric radius 
defined through the Sachs form factor Ge- The quantity Oo calculated in a non-relativistic 
quark model without taking into account relativistic corrections leads to the same magnitude 
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of cto for proton and neutron. Since Aa is equal to zero for the neutron but gives a significant 
contribution to a for the proton, one can naively expect that 

«p — ttn = = 3.8 ± 0.1 . (6) 

Comparing (§) with we would come to the conclusion that the experimental and predicted 
differences between the electric polarizabilities of proton and neutron deviate from each 
other by 2.3(T. In the following it will be shown that this discrepancy may be connected 
with hitherto unknown relativistic corrections to Oo which are of the same order as Aa. 
The relativistic corrections to Oo in (^) come from corrections to the wave functions, to the 
energies of the ground and excited states, and a correction to D. The relativistic correction 
to the electric dipole operator is connected with the appropriate relativistic definition of the 
center-of-mass coordinate [ll3| , |14| , |l5l . The neglect of this correction leads to an incomplete 



expression for a, and, as it has been shown in |15], the missing piece can be very essential. 



In the present paper we show that the difference 5a = (ao)p— (ao)n is due to the correction 
to the electric dipole operator only. By carrying out model calculations we show that 6a can 
strongly compensate Aa, thus removing the contradiction between (j^) and (|^) mentioned 
above. 



2 Correction to the center-of-mass vector 

Let us consider a system of spin- 1/2 particles having the masses and spin operators Sj. 
The Hamiltonian of the system H , accounting for the first relativistic correction, is the sum 
of the non-relativistic Hamiltonian Hnr and the Breit Hamiltonian Hb, where 

i * i>j 

Here = — Tj and Vij{r) are some potentials. The explicit form of the Breit Hamiltonian 
Hb depends on the Lorentz structure of interaction, namely, whether it is Lorentz-scalar 



or Lorentz- vector (see |jT6[ and references therein). As in most models, we assume that 
the potential V{r) contains a short-range Coulomb-type term Vc{r) due to single-gluon 
exchange corresponding to vector interaction, and a long-range confining potential Vconfi^)- 
The Lorentz structure of the latter is not well known. We show that it makes no difference 
in the form of 6Tir m ■ 



The center-of-mass vector Rem. should satisfy the following relations 

P P / if \ 

[R,^ , Hnr + Hb] = i— ~ «TT 1 - ^ . (8) 

L cm. , nr Bi M + H M \ M J ^ ^ 

where M and P are the total mass and the total momentum of the system, respectively. 
Using the explicit form of the Breit Hamiltonian it is easy to check that in both cases, i.e. 
Lorentz-scalar and Lorentz-vector interactions, the vector Rem. for the system of particles 
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has the form 



i 



2M ^ V ^ 



R, 



X s, 



Pi 



rrii 



(9) 



where M = ^mj, and where the notation {a,b} = ab + ba is used. Thus, the relativistic 
correction SHq is expressed via the non-relativistic potential of the interaction between the 
constituents. 

Let us consider proton and neutron. In this case all potentials are equal to each other, 
i.e. Vij = V. Let the vectors ri and r2 correspond to two w-quarks in the proton and two 
d-quarks in the neutron, and correspond to the d-quark in the proton and the n-quark in 
the neutron, and mi to the mass of the identical quarks {u quarks in the proton and d-quarks 
in the neutron), and to the mass of the third quark. We pass to the Jacobi variables: 

11 1 
P = ;^(ri-r2) , A = -^(ri + rs - 2r3) , R = — (miri + mirs + mgra) , (10) 

where M = 2mi + m^. Then the momentum operators are 

1 ^/6 
Pp = ;^(Pi - P2) , Pa = ^("^3Pi +"^3P2 -2miP3) , P = Pi + P2 + Ps • (H) 

The non-relativistic Hamiltonian H„r in terms of the Jacobi variables has the form 



+ 



Pp 



Px 



2M 2mi 2mx 



+ V{V2p) + V{V20 + V{V2r]) , 



(12) 



where mx = Smim-^/M, ^ = (p + \/3A)/2, and 77 = (— p + -\/3A)/2. In terms of Jacobi 
variables the spin-independent correction to the center-of-mass vector is 



SRn 



2M2 



({p,P-p,} + {A,P-p,}) 



2V6miM 



1 V6 
Ip , Pa • Pp} 



M2 



>^ , — Pp + Pa 



^^VemsiVM + 



2M2 



p A/6(m3 - 2mi)A 
y/2 2M 



2M 



+ 



p y6(m3 - 2mi)A 
\/2 2M 



2M 



(13) 



and the spin-dependent correction 



P X S Pp X (si - S2) _^ Pa_x [m3(si + S2) - 2miS3] 



2M2 2^/2mlM 



2-\/6mim3M 



(14) 
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When mi = 
essentially simpler: 



m , then m\ = m and the expressions for SHq and SUg become 



and 



SRn 



[{p,P-p,} + {A,P.p,}] 



10 o ILf 1 ^ sr'p 



1 



1 



m 



{p, PA-Pp} + -{A,pJ-p^} 
2\V{V2p) + (VSp - X) V{V20 + (-VSp - X) V{V2r]) 



SR. 



1 1 

3P X S + -j=Pp X (si 



1 

S2) + ^PA X (Si + S2 

V6 



2s, 



(15) 



(16) 



Let us consider the substitutions 

1 Vs. 

p ^ -p A 

^ 2^^ 2 



1 



P/^ 



V3 



Pa 



Pa 



^PP + 



2PA^ 



(17) 



If mi = ms , then the non-relativistic Hamiltonian ([I4) is invariant under this substitution 
and 5Ro(P = 0) ^ — (5Ro(P = 0). In addition, the Hamiltonian (0) is invariant under the 
transformation p —>■ —p and A — > —A. The properties of the operators with respect to the 
transformations (IT^) are useful for the selection rules for the matrix elements. 



3 Charge radii 

Let us consider the charge radii of the nucleon r"^ = (Yl ~ ^c.m.Y ), where e, is equal 
to 2/3 for the w-quark and —1/3 for the (i-quark: 

= ^(0| [4(ri-Re.™.)2-(r3-Re.™.)'] |0), 
{rl)n = -^(0| [(ri-R,.,„.)2_(r3_R_)2] |o). (18) 

Here we used the symmetry of the wave function with respect to the permutation ri ^ r2. 
The values of {r%)p and {r'^)n were measured to be {r'^)p = 0.74 fm^ and {r%)n = —0.119 ± 
0.003 fm^ 113. Thus, (0|(ri - Rc.^.)2|o) > (0|(r3 - Rc.,„.)2|o). Due to the symmetry of i/„^ 



with respect to the permutation ri ^ for the case of equal quark masses, it is evident 
that {rE)n = in the non-relativistic approximation. The relativistic correction to {r'^)n 
comes from the correction to the center-of-mass vector and from the correction to the wave 
function: 

S{rl)n = ^(0| {ri - ra , <5Ro + 5RJ |0) 
-^(0| [(n - R)2 - (r3 - R)2] GoHb\0) 

-1{0\HbGo [(ri - Rf - (rs - R)^] |0) , (19) 
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where Gq is the non-relativistic reduced Green function 

Go = [eo-Hnr + tO]-\l-\0){0\), (20) 
and £0 is the ground state binding energy in the non-relativistic approximation. 

Using the properties of 5Ro and 6Hs with respect to the transformations (p!7|) , it is 



easy to show that the contribution to ^{r'^)n ( P^ from the correction to the center-of-mass 
vector, as well as the contribution of the spin-independent part of Hb, vanish. Due to parity 
conservation the contribution of the part of Hb being linear in spin is also zero. The spin-spin 
part H^^^ of the Breit Hamiltonian gives a non-zero contribution to 5(r|;)„. This operator 
exists only for the Lorentz- vector part Vu{r) of the potential: 

- E 7^;^.^^"^'''^^ - K(r..)M,][3(s, ■ f,,)(s, ■ f,,) - s, ■ s,] , (21) 

where = Tij/rij. Substituting Vj;{r) = —2as/'ir and averaging over the spin part of the 
neutron wave function, we obtain 

H^''^ = [Kp) - m) - nv)] . (22) 

Using Jacobi variables, we can represent 5(r|;)„ as 

S{rl)n = -Imp' - \']GoH^''^\0) - i(0|/7(^^)Go[/ - A^]|0) . (23) 



In order to calculate the matrix element (^) we follow the prescription of \TE\ and set 
V{r) = Kr'^/2 in the non-relativistic Hamiltonian Hnr- Then we obtain 

H ^EL + ^ + A + ^JSil±2^ (24) 
2M 2m 2m 2 ' ^ ^ 

Thus, we have two independent oscillators with equal frequencies ujp = ux = ^^/3K/m. 
Taking into account that (p^ — A^)|0) is an eigenfunction of Hnr ( p4D with the excitation 
energy E — Eq = 2ujp, we have 

= T^(0|[p^ - \']H^-^\0) = - . (25) 

Thus, we obtain the negative value for (r|;)„. The magnitude of this quantity can also be 
made in agreement with the experimental value by taking the appropriate parameters. 



4 Electric polarizability 

Using the Jacobi variables we obtain the operator of the internal dipole moment D = 
e ^ ej(rj — Rem.) for proton and neutron: 
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If we neglect SHc.m. in Dp, we immediately obtain (a;o)p = (ao)n since ao is quadratic in D 
(see (^)). Therefore, the difference between (ao)p and (ao)n arises only due to the correction 
SHc.m. in Dp. The contribution to this difference being hnear in 6Hc.m. reads: 

6a = («o)p - (ao)n = f 3 j [(0|AG(5R,.„.|0) + (0|5R,.„.GA|0)] , (27) 

where G is the reduced Green function accounting for the first relativistic correction: 

G=[Eo- H„r -Hb + tO]-\l - |0)(0|) . (28) 

If we replace G by its non-relativistic limit Gq and use the symmetry with respect to the 
transformations (|l3), we obtain zero as the result for 6a. Therefore, it is necessary to take 
into account the corrections to the wave function and to the Green function due to the 
spin-dependent part of Hb- In addition, it is necessary to account for the term quadratic in 
6T{.c.m. in (^), and the second-order relativistic correction to Rem.- The calculation of the 
latter is a very complicated problem. Since we are only going to demonstrate the possible 
cancellation between 6a and Aa (see (^)), we may simplify our problem and assume, in a 
spirit of the diquark model that the relativistic effects reduce to the small difference 
between the masses mi and m^. As a result, we use the non-relativistic Hamiltonian 

Hnr = + + + (29) 

2M 2mi 2mA 2 

and the expression ([T3| ) for ^Rq. Remind that M = 2mi + and m^ = Smim^/M. Thus, 
the corresponding frequencies are Up = ^J'iKjmx and uj\ = ^J^Kjmx- We can fix the 
parameters of the model from the experimental values of [f^Blp ^"^^ ('^l;)"- Using (p!0D , (|18D 
and (P^D it is easy to find that 



E}v 



EJn 



rriiujp 
1 



2miu; 



m^Mojx 
(2mi - mz)ujp 



miojp 



p L 



5(1-2;) 
Grriiujp 



(30) 



where x = mi/m^. Substituting the experimental values of charge radii into (^), we obtain 
X = 0.8 , and rriiUp = 6 ■ lO^MeV^. Using the conventional value mi = 330 MeV, we have 
^p/fT^i = 0.6. As a result of simple calculations we have in our model: 

e- 79e2(l-x) 536^ 

Aar, = 7, — , 6a = ^ 1 5- . (31) 

^ 9mjujp WSmjujp 1296m? ^ ^ 

The second term in 6a comes from the correction being quadratic in 6Iic.m. and numerically 
is much smaller than the first one. For the parameters of our model we obtain 

^^±^ = -0.1, (32) 
Aap 

which leads to the approximate equality of the predicted proton and neutron electric polar- 
izabilities in agreement with the experimental data. 
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Conclusion 



Thus we have demonstrated that the equahty of the proton and neutron electric polarizabihty 
can be explained in the frame of constituent quark model. Therefore, it is not necessary to 
deal with such ad hoc contributions (for constituent quark model) as mesonic currents. 
Though it has been shown in that there are many sources of the relativistic corrections 
of the same order as Aa, due to the isospin invariance it is crucial to take into account 
the correction to the center-of-mass vector. If this correction is neglected, the difference 
(«o)p — {cto)n vanishes identically. The estimate for the relativistic corrections made above 
with the use of a toy model serves only as an illustration of the possibility to make the 
difference {a)p — (q;)„ consistent with the experimental data at reasonable values of the 
parameters. 
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